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SUMMARY 

/ 

In  the  past  6  months,  research  has  concentrated  on  the  development  and 
application  of  a  surface-wave  magnitude  scale  based  on  the  amplitude  of  nigher- 
mode  Rayleigh  waves.  A  relatively  high  ratio  of  higher  mode  to  funda¬ 
mental  mode  energy  can  be  produced  when  an  unusual  combination  of  source  depth 
and  focal  mechanism  leads  to  poor  excitation  of  the  fundamental  mode  Rayleigh 
wave,  in  which  case,  the  amplitude  of  the  higher  mode  yields  a  more  reliable 
estimate  of  the  size  of  the  event  than  the  amplitude  of  the  fundamental  mode. 
Several  anomalous  events  in  central  Asia,  which  are  characterized  by  unusually 
low  Ms :mb  ratios  and  hence  could  be  suspected  to  be  nu-lear  explosions,  are 
reclassified  as  earthquakes  when  the  new,  higher-mode  scale*  Mg»  is  employed. 

The  focal  depth  and  mechanism  of  these  anomalous  events  apparently  causes  poor 
excitation  of  the  fundamental  mode  Rayleigh  wave  without  significantly  affecti  g 
the  amplitude  of  the  higher  mode.  Elimination  of  these  source  related  factors 
from  traditional  Ms:n>b  discriminants  is  an  essential  tool  for  reliable 
discrimination  between  earthquakes  and  nuclear  underground  tests. 
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I.  STATION  MAINTENANCE 

The  high  gain  station  at  Ogdensburg,  New  Jersey  has  been  maintained  in 
constant  operation.  Approximately  90  photographic  seismograms  are  acquired  each 
week  in  addition  to  a  magnetic  tape  once  per  three  weeks.  These  data  are 
forwarded  to  the  Albuquerque  Seismological  Center  and  a  copy  of  the  tape  is 
stored  in  the  Lamont  Data  File.  The  displacement  outputs  of  the  Geotech 
long-period  seismometers  are  recorded  on  2"  per  hour  analogue  chart  recorders 
as  are  the  outputs  from  the  four  quartz-tube  strainmetcrs .  We  are  in  the  process 
of  transferring  the  seven  chart  records  to  a  multichannel  3"  per  hour  recorder 
which  presently  records  data  from  two  new  wire  strainmeters . 

The  routine  operation  of  the  observatory  has  required  two  fu  1  time 
technicians  in  the  past.  From  October  1975  we  have  operated  the  station  with 
one  technician  (M.  Connor)  and  have  discontinued  recording  seismic  data  from 
levels  others  than  the  -560  m  level  of  the  mine  in  which  the  observatory  is 

located • 

The  digital  recording  of  the  long-period  data  has  been  interrupted  on  a 
number  of  occasions  by  electronic  malfunction  of  the  Astrodata  recorder. 

This  is  now  maintained  by  the  Albuquerque  Seismological  Center  who  have 
arranged  on  site  repairs  when  required. 

II.  RESULTS  OF  THE  DATA  ANALYSIS 

A.  Higher  modes  as  an  aid  to  seismic  discrimination, 

A  paper  describing  the  application  of  a  surface  wave  magnitude  scale  based 
on  the  amplitude  of  higher-mode  Rayleigh  waves  to  the  discrimination  of 
anomalous  earthquakes  from  underground  nuclear  explosions  has  been  completed 
and  accepted  for  publication  in  the  Bulletin  of  the  Seismological  Society. 

higher-mode  scale,  M^,  several 


of  America  (Forsyth,  1976).  By  using  the  new 


3 


anomalous  events  in  central  Asia  reported  by  other  investigators  were  shown  to 
be  normal  earthquakes.  The  conclusions  based  on  the  higher  modes  were  confirmed 
in  a  number  of  cases  by  the  observation  of  long-period  (AO  to  60  s)  fundamental 
mode  signals  at  the  high-gain  long-period  station  at  Chiang  Mai.  Thailand.  A 
relatively  high  ratio  of  higher  mode  to  fundamental  mode  energy  can  be  produced 
when  an  unusual  combination  of  source  depth  and  focal  mechanism  leads  to  poor 
excitation  of  the  fundamental  mode  Fayleigh  wave,  in  which  cat-,  the  amplitude 
of  the  higher  mode  yields  a  more  reliable  estimate  of  the  sire  of  the  event  than 
the  amplitude  of  the  fundamental  mode.  Several  anomalous  events  in  central 
Asia,  which  are  characterized  by  unusually  low  Ms:mb  ratios  and  hence  could 
be  suspected  to  be  nuclear  explosion...  are  reclassified  as  earthquakes  when  the 
new.  higher-mode  scale.  M*.  is  employed.  The  focal  depth  and  mechanism  of  these 
anomalous  events  apparently  causes  poor  excitation  of  the  fundamental  mode  Rayleigh 
wave  without  significantly  affecting  the  amplitude  of  the  higher  mode.  Elimination 
of  these  source  related  factors  from  traditional  M8=mb  discriminants  is  an  es¬ 
sential  tool  for  reliable  discrimination  between  earthquakes  and  nuclear  under- 

ground  tests* 

There  are  three  primary  cases  in  which  a  relatively  high  ratio  of  hiR 
mode  to  fundamental  mode  energy  will  be  generated  at  the  source.  1)  The 
earthquake  may  be  below  the  depth  at  which  the  fundamental  mode  is  excited 
efficiently,  yet  still  be  sufficiently  shallow  to  excite  higher  modes  within 
the  crustal  wave  guide.  2)  For  certain  focal  mechanisms,  fundamental  mode 
Rayleigh  waves  will  be  excited  only  weakly  if  the  earthquake  occurs  at  a  depth 
corresponding  to  a  zero  crossing  in  displacement  or  stress  versus  depth  for 
the  wave.  Zero  crossings  for  the  higher  modes  are  generally  not  at  the  same 
depths,  so  if  the  amplitude  of  the  fundamental  mode  is  reduced  by  this  effect. 
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the  higher  nodes  nay  be  undiminished.  3)  Except  for  sources  very  near  the 
surface,  the  azimuthal  radiation  pattern  for  higher  nodes  may  differ  from  the 
fundamental  nodes.  At  azimuths  corresponding  to  a  node  in  the  fundamental  mode 
radiation  pattern  there  may  appear  to  be  a  relatively  large  amount  of  higher 
mode  energy.  Although  each  of  these  three  mechanisms  for  producing  a  large 
ratio  of  higher  mode  to  fundamental  mode  energy  is  frequency  dependent,  they 
can  affect  a  range  of  frequencies  broad  enough  to  significantly  change  the 
appearance  of  the  seismogram.  All  three  mechanisms  involve  reducing  the  ampli¬ 
tude  of  the  fundamental  mode  rather  than  increasing  the  absolute  amplitude  of 
the  higher  modes.  Thus,  when  the  higher /fundamental  ratio  is  large  and  the 
surface  wave  magnitude  based  on  fundamental  mode  Rayleigh  wave  amplitudes 
is  anomalously  low  compared  to  the  body  wave  magnitude,  a  more  reliable  esti¬ 
mate  of  seismic  moment  or  surface  wave  magnitude  may  be  obtained  from  the 


amplitudes  of  the  higher  modes. 

As  one  of  the  primary  applications  of  an  improved  magnitude  scale  will  be 
to  aid  in  discrimination  between  earthquakes  and  underground  explosions,  the 
higher  mode  scale  (M*)  is  designed  to  be  equivalent  to  the  fundamental  mode 
scale  (Mg)  for  a  near-surface  focus  event.  The  scale  now  commonly  accepted  for 
Asian  events  (liar shall  and  Basham,  1972;  Conference  of  the  Committee  on 

Disarmament,  1972)  is  of  the  form 

Mg  °  log  A  +  B*  (A)  +  P(T)  ^  ^ 

where  A  is  the  maximum  amplitude  (nm)  of  the  vetticel  component  of  ground  motion 
in  the  Rayleigh  wave  train,  T  is  the  period  of  the  wave  of  maximum  amplitude, 

P(T)  is  the  period  correction  which  depends  on.the  propagation  path  and  B'(A)  is 
a  distance  normalizing  term  that  corrects  for  the  effects  of  geometric  spreading, 
scattering  and  absorption  of  the  propagating  surface  wave.  P(T)  is  assigned  value 
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0.0  at  the  arbitrary  reference  period  of  20s. 

To  equalize  Mfl  and  M*,  it  is  necessary  to  know  the  relative  amplitudes  of  the 

two  modes.  At  the  surface,  the  excitation  function  reduces  to  a  simple  form 
with  only  those  terms  proportional  to  the  horizontal  displacement  remaining. 

The  amplitude  of  the  fundamental  mode  wave  of  20s  period  should  be  approximately 
7.5  times  the  amplitude  of  the  first  higher  mode  wave  of  10s  period  for  any 
very-near-surface  focus  event.  The  formula  for  Mj  should  therefore  be: 

Mg  -  log  A  4*  B *  (A)  +  log  T/T0  +  log  7.5 

where  T  is  the  period  of  the  maximum  amplitude  A,  B'(A)  is  distance  correction 
from  Marshall  and  Basham  (1972),  and  TQ  is  the  reference  period  of  10s.  This 
formula  should  be  modified  if  used  in  areas  other  than  central  Asia,  because 
B' (A)  and  the  mode  correction  factor,  log  7.5,  may  vary  from  region  to  region. 

The  ratio  of  surface-to  body-wave  magnitude,  Ms:mb,  provides  one  of  the 
principal  means  of  distinguishing  earthquakes  from  underground  explosions  (eg. 
Evernden,  1969;  Liebermann  and  Pomeroy,  1969).  Explosions  tend  to  be  less 
efficient  generators  of  surface  waves,  with  Mg  less  than  [tnb  -  1.5].  In  con¬ 
trast,  M  for  most  shallow  earthquakes  is  greater  than  or  equal  to  [mb  -  1.0] 
(Conference  Committee  on  Disarmament,  1972).  However,  there  are  a  few 
"anomalous”  earthquakes  which  have  explosion-like  Mg:mb  ratios,  or  fall  between 
the  earthquake  and  explosion  population,  i.e.  [mb  -  1.0]  >  Mg  >.  [mb  -  1.5],  and 
thus  are  statistically  indistinguishable  from  explosions.  Examples  of  anomalous 
earthquakes  in  central  Asia  have  been  presented  by  Landers  (1972),  the  Conference 
of  the  Committee  on  Disarmament  (1972),  Douglas  et  al.  (1974)  and  Nuttll  and 


Kim  (1975). 
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The  long-period,  vertical  component  of  motion  for  two  anomalous  earth¬ 
quakes  and  an  underground  explosion  are  shown  in  figure  2  as  recorded  at  the 
WWSSN  station  in  Kabul,  Afghanistan  (KBL,  Fig.  1).  Locations  and  magnitudes 
are  given  in  Table  1.  The  signal  from  event  A,  an  earthquake  in  Kazakh,  shown 
at  the  top  of  the  figure,  happens  to  be  superimposed  on  the  Rayleigh  wave  train 
from  a  much  larger  earthquake  in  the  Tonga  region.  Although  the  high-frequency 
Rayleigh  wave  from  the  Kazakh  event  is  clearly  visible,  it  is  impossible  to 
obtain  a  reliable  estimate  of  the  20s  signal  level  of  the  smaller  earthquake. 
Event  A  can  be  seen  at  several  other  stations  in  central  Asia,  but  in  each 
case,  the  surface  waves  from  Tonga  interfere.  All  of  the  anomalous  events  are 
small,  making  it  difficult  to  measure  the  spectrum  over  a  period  range  broad 
enough  to  avoid  frequency-dependent  source  effects,  even  when  there  is  no 
interfering  signal.  For  event  A,  Ms  has  to  be  based  on  signals  of  period  7  to  8s 
yielding  an  estimated  magnitude  of  2.9,  well  within  the  explosion  population. 
However,  as  Landers  (1972)  pointed  out,  the  amplitude  of  the  higher  mode  arrival 
at  about  3.4  km/s  is  equal  to  or  greater  than  the  amplitude  of  the  fundamental 
mode.  The  revised  estimate  of  surface  wave  magnitude,  =  3.9,  places  the 
event  in  the  earthquake  population  (Table  1).  Another  similar  event  located 
nearby  (C  in  Table  1)  is  also  removed  from  the  anomalous  category  by  using  Mg. 
The  number  of  observations  on  which  M*}  is  based  is  listed  in  the  last  column 

of  Table  1. 

Unlike  events  A  and  C,  for  which  no  long-period  surface  waves  could  be 
detected,  event  B  generated  20s  Rayleigh  waves  which  could  easily  be  measured 
at  KBL  (2nd  line  in  fig.  2).  The  M8*.mb  relation  for  event  B  falls  between  the 
earthquake  and  explosion  populations,  but  again,  the  amplitude  of  the  higher 
mode  restores  it  to  the  earthquake  group.  m£  also  allows  two  of  the  anomalous 
earthquakes  listed  by  Nuttli  and  Kim  (1975)  to  be  classified  as  normal  earth- 
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qaakcs  (events  E  and  F  in  Table  2).  Event  F  is  unusual  in  that  40-to  60s  petiod 
Rayleigh  waves  ate  clearly  visible  at  the  »SSH  station  HM  and  the  hi8h-gain. 
long-period  station  at  Chiang  Mai,  Thailand  (C1IG) ,  even  though  very  little 
fundamental  mode  Rayleigh  or  hove  wave  signal  with  period  less  than  30s  can  be 
seen  at  any  station.  The  higher  mode  from  this  event  travels  with  a  group  ve¬ 
locity  of  about  3.8  km/s,  in  contrast  to  the  crustal  velocities  observed  for 
other  events.  The  mantle-like  velocity  of  the  higher  mode  and  the  long-period 
character  of  the  fundamental  mode  suggests  that  the  earthquake  occurred  deep 
within  the  Himalayan  crust,  perhaps  at  a  depth  of  70  km  or  more.  Ms  based  on 
the  60s  period  wave  is  about  6.0,  in  reasonably  good  agreement  with  mJ.  Event 
D  occurred  in  nearly  the  same  location  as  event  F,  but  was  too  small  to  generate 

an  observable  higher  mode.  The  fundamental  mode  could,  however,  be  detected 

a  4 <-  «=  recorded  with  the  same  waveform  at 

at  the  CHG  high-gain  station  where  it  was  recorded 

about  half  the  amplitude  of  event  F.  1  suggest  that  event  D  also  occurred  deep 
within  the  Himalayan  crust  and.  by  analogy  with  event  F,  should  be  assigned  an 
Hh  of  about  3.8.  Other  examples  of  higher  modes  from  anomalous  earthquakes  in 
L  Tibetan  area  are  presented  by  Tatham  at  al.  (1975)  and  the  effect  of  using 
instead  of  Mg  is  shown  in  fig.  3. 

If  the  higher  mode  magnitude  scale  is  correct,  M*  should  be  equal  to  M, 
for  very  shallow  earthquakes  or  underground  explosions.  Although  there  are  as 
yet  no  known  explosions  within  an  area  with  a  70  km  thick  crust,  an  attempt  was 
made  to  test  the  validity  of  the  scale  using  signals  from  the  Soviet  nuclear 
test  site  in  eastern  Kazakh  (fig.  1).  One  effect  of  a  thinner  crust  is  to  reduce 
the  period  of  the  higher  modes  trapped  in  the  crustal  wave  guide.  As  can  be 
seen  in  the  bottom  trace  in  fig.  2,  the  maximum  period  of  the  higher  mode 
generated  by  the  explosion  is  only  atout  5s.  +  for  three  explosions  at  the 

site,  events  0,  H  and  1  in  Table  2,  is  0.05  to  0.15  units  larger  than  K.. 


This 
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discrepancy  is  nor  large  and  nay  be  dne  to  the  dlficreaces  in  crustal  structure 
between  the  test  site  and  the  locations  of  the  anomalous  earthquahes  for  which 

the  scale  was  designed. 

B.  Ultra-long  period  studies 

The  Jong-pcriod  analogue  records  from  OGD  strainmetets,  gravitetets  and 
tiltmeters  accumulated  at  the  seismic  observatory  since  1,60  are  being  searched, 
cataloged  and  reorganised.  The  data  is  in  the  formof  strip  charts  and  a 
machine  has  been  devised  to  enable  the  rapid  conversion  of  the  analogue  trace  to 
a  digital  card  format.  The  strain  records  are  of  Interest  in  that,  although 
they  are  poorly  calibrated  in  places  and  possess  uncertain  linearity,  they 
represent  an  unusually  long  data  set.  We  intend  to  search  the  data  for  temporal 
changes  in  relative  tidal  admittance  and  to  examine  the  ultra-long-period  part 
of  the  seismic  spectrum.  The  re-organiration  revealed  that  the  guartt  strain- 
meters  have  operated  for  ,5%  of  the  time  since  their  installation,  the  missing 
5*  being  partly  by  the  loss  of  individual  rolls  of  data  which  have  been  borrowed 
and  not  returned.  Except  for  periods  of  maintenance  and  occasional  electronic 
malfunction  most  of  the  existing  data  are  usable.  The  strain  data  consist  of 
approximately  1100  rolls  of  paper  chart  with  a  single  ink  trace  written  on 
each.  During  the  reorganisation  each  roll  was  examined  for  time  math  consistency. 

calibration  marks,  manual  adjustments  and  for  data  quality. 

Some  of  the  strainmeter  data  has  already  been  digitized.  A  plot  of  4 
years  of  digitized  data  from  one  strainmeter  is  shown  in  fig.  4  and  its  stability 
in  comparison  with  other  high-quality  instruments  is  shown  in  fig.  5.  The 
record  has  been  smoothed  with  a  48  hour  Gaussian  Filter  to  remove  tides  «, 

t  \  Thp  record  reveals  two  important  features 
tidal  amplitude  shown  to  scale). 

4  daft,  There  is  an  absence  of  annual  thermal  signals 
of  the  Ogdensburg  strain  data.  Ther 

,  ifi-7/„onr  We  believe  that  at  least  some  of  the 
and  the  "drift"  rate  is  less  than  10  /year.  We  ben 
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observed  data  Is  affected  by  instrumental  instability  since  increased  otrai 
rates  follow  instrumental  modifications.  It  appears  that  the  record  may  be 
of  sufficiently  low  instrumental  drift  to  comment  eventually  on  changes  of 

secular  strain. 

XXI.  improvements  to  long-period  recording  systems 

The  Ogdensburg  strainmeters  have  been  recorded  since  1960  on  paper  chart 
analogue  recorders  (at  2  inches  per  hour  without  filtering)  and  as  seismograms 
(at  *  inch  per  minute)  after  passing  through  a  band-pass  filter.  The  dynamic^ 
range  of  the  band-pass  filter  is  inadequate  for  useful  studies  of  the  Earth's 
free  oscillations  and  the  characteristics  of  the  filter  are  insufficient  to 
suppress  20  second  surface  waves  relative  to  waves  with  periods  between  1 
minute  and  1  hour.  The  circuit  diagram  of  the  filter  is  shown  in  fig.  6. 

We  intend  to  improve  the  recording  system  to  enhance  the  dynamic  range 
and  accuracy  of  free  oscillation  data  both  from  the  strainmeters  and  the 
Broad-band  seismometer.  Ultimately  we  shall  record  four  components  of  strain 
(azimuths  29.5® ,  69.5®,  132®  and  vertical)  and  three  components  of  acceleration 
(the  displacement  outputs  of  Geotech  N/S,  E/W  and  vertical  seismometers)  on 
a  digital  recorder  sampling  at  5  second  intervals.  As  an  immediate  step 
towards  this  goal  we  have  introduced  Improved  filters  to  replace  the  previous 
band-pass  filters.  The  circuit  diagram  and  output  charaC  eristics  of  the  new 

,  ,  . v  u  E  Farrell,  CIRES  Boulder)  are  shov.i  in  figs.  6  and  7. 

filter  (designed  by  w.E.  tarre-n, 

The  filters  have  been  modified  to  interface  with  present  recording  equipment 
and  have  been  provided  w'th  differential  high-impedance  Inputs.  Seven  filters 
are  housed  in  a  hermetically-sealed  enclosure  which  now  acts  as  a  distribution 
point  for  all  the  long-period  strain  and  acceleration  data. 

Bering  the  coming  year  we  will  upgrade  the  displacement  transducers  on  all 

At  present  these  are  all  tuned-resonant  circuit  capacitative 


seven  instruments. 
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systems  which  have  poor  stability  and  are  difficult  to  calibrate.  We  intend 
to  replace  them  with  ratio-type  transformer  capacitive  transducer  electronics 
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Fig.  1 

Fig.  2 


Fig.  3 

Fig.  A 

Fig.  5 


FIGURE  CAPTIONS 

Index  map  of  source  locations  and  WWSSN  stations.  Events  1-3 
are  listed  in  Table  1  and  events  A-I  are  listed  in  Table  2. 

Tracings  of  long-period,  vertical  component  of  records  at  the 
WWSSN  station  in  Kabul  (KBL) .  Traces  have  not  been  normalized 
and  relative  amplitudes  are  correct.  Arrows  indicate  group  velocities 
on  each  trace.  Top  line  is  from  anomalous  event  A,  11°  away; 
middle  line  is  from  anomalous  event  B,  6°  away;  and  bottom  line  is 
from  underground  explosion  G,  17°  away.  Strong  arrival  at  about 
3.0  km/s  on  explosion  record  is  short-period  branch  of  fundamental 

mode. 

Mg-rn^,  diagram  summarizing  analysis  of  Tibetan  anomalous  events. 

Arrows  and  diamonds  show  increase  in  Mg  values  when  observations  of 
higher-mode  Rayleigh  waves  are  included,  with  (+)  indicating 
original  Ms  value.  Arrows  and  open  circles  show  increase  in  Mg 
values  when  observations  of  Love  waves  are  included.  All  events 
are  located  in  eastern  Tibet  and  are  described  in  detail  in  Tatham 

et  al.  (1975). 

Ogdensburg  strain  data  with  tides  removed  by  A8-hour  Gaussian 
filtering.  M2  tide  amplitude  shown. 

Ogdensburg  secular  strain  compared  with  strain  rates  observed  by 
other  investigators.  The  Pinion  Flats  data  when  corrected  for 
tilting  of  end  piers  exhibits  considerably  less  annual  strain.  The 
data  is  presented  in  "year s-since- installation".  Strain  in  units  of 

10'7. 

Circuit  diagrams  for  old  (bottom)  and  new  (top)  band-pass  filters 
for  strainmeter  data.  New  filters  are  more  stable,  with  flatter 


Fig.  6 
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response  and  greater  dynamic  range. 

Fig.  7  Amplitude  and  frequency  response  of  outputs  from  new  strainmeter 
filter. 


Figure  2. 
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